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Rua de Alexandre Herculano, 244 — PORTO 
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DELEGAÇÃO 
Rua Rosa Araújo, 35 — LISBOA 
TEL.: 731181/82/83 81/85/86 - 733333 - 18491 


ELECTRICIDADE DO LINDOSO 
E DAS CENTRAIS DO FREIXO 
E DA CACHOFARRA 


a 
A UNIÃO ELÉCTRICA PORTUGUESA distribui, 


pela mais extensa rede de alta tensão em Por- 
tugal (cerca de 2500 km), energia eléctrica em 
70 concelhos, nos distritos de VIANA DO CAS- 
TELO, BRAGA, PORTO, AVEIRO, COIMBRA, 
VISEU, LEIRIA, SETÚBAL e ÉVORA. 


Oficina e Laboratórios 
do 


INSTITUTO SUPERIOR TÉCNICO 


As oficinas pedagógicas do Instituto 
Superior Técnico, de CARPINTARIA 
DE MOLDES, de INSTRUMENTOS 
DE PRECISÃO de ELECTROTE- 
CNICA, fornecem todo o género de 
material escolar de demonstração 
para o ensino técnico. Nos laboratórios 
de QUÍMICA - ANALÍTICA, FÍSICA 
INDUSTRIAL E DE MINERALOGIA 
executam-se análises para o público. 


Para quaisquer informações, dirigir-se 
ao secretário da comissão executiva 
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LAVA-LOUÇAS EM AÇO INOXIDAVEL 


Modermze os lares de Portugal, usando mas suas com 
nhas Lavalouça em AÇO INOXIDAVEL. Concepção 
“ execução da 


FABRICA PORTUGAL 


elegantes e higiênicos - duração ilimitada - revestidos de isolante 
especial - fornecidos com sifão de segurança - não partem a louça 
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PIMENTEL & 
CASQUILHO, L.” 


Quase 40 ANos 
AO SERVIÇO DA 
ENGENHARIA PORTUGUESA 
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PRECISÃO, CIENTÍFICOS 
E INDUSTRIAIS 
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Bombas centrífugas «PACIFIC» 


Construídas por 


Batignolles = Chatillon 


Especialmente destinadas para: 


— Petróleos 

— Indústrias quimicas 
— Construções navais 
— Centrais eléctricas 


— Alimentação de caldeiras 


— Etc., etc. 


Sob as condições de funcionamento 
seguintes: 


— Débitos: 1 m*/h a 10.000 m*/h 
— Pressões: altas e baixas 

— Temperaturas: — 180º C + 600ºC 
— Líquidos: corrosivos ou não 


Para estudos e orçamentos, consultar o representante em Portugal: 
CIDEX — Comércio Internacional 
e Representações Industriais, Lda. 
Av. Casal Ribeiro, 46-2.º/Telef. 7353 07/LISBOA1 


[X 
SONDAGENS/E: 


55 
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MOTORES PARA TODAS AS APLICAÇÕES INDUSTRIAIS 
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OUTROS MODELOS PARA ENTREGA MUITO RÁPIDA 


MOTORES A 4 TEMPOS - ARREFECIMENTO POR ÁGUA 
ARRANQUE A AR COMPRIMIDO OU ELÉCTRICO 


REPRESENTANTES GC. SANTOS, S.A. R. L. - pivisão MARÍTIMA E INDUSTRIAL 


29, AV. DA LIBERDADE, 41 -LISBOA 
160, R. DE S.tTA CATARINA. 168-PORTO 
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EM CADA CHAMADA PRESTAMOS UM SERVIÇO 


O seu telefone, servo obediente, aprende novas maneiras de o servir. Muitas das vantagens já 
hoje existentes nas centrais telefónicas internas serão num futuro próximo introduzidas nas 
redes públicas. Assim, sem ter que levantar o auscultador, poderá marcar o número preten- 
dido. Se por coincidência este estiver interrompido, o seu telefone avisá-lo-á do momento em 
que a linha está livre evitando, assim, grandes perdas de tempo e aborrecimentos. Por outro 
lado a telefonista pode fazer uma chamada e efectuar imediatamente a sua transferência sem 
ter que preocupar-se mais. Esse sistema actua como que se gravasse o que lhe foi dito: será, 
pois, um fiel servidor. Estas e outras extraordinárias realizações já são incluídas nos equipa- 
mentos da ITT. 

Centrais de tipo convencional que ocupam salas darão lugar a equipamentos com dimensões 
de gavetas. 

Estes projectos tipificam o carácter pioneiro da ITT em matéria de comunicações. Já em 1930 
a ITT havia instalado um sistema de chamada directa entre assinantes a longa distância. Traba- 
lhos desta natureza, bem como muitos outros aplicados à TV, radar, micro-ondas, comunicações 
marítimas e espaciais, permitiram a ITT (International Telephone and Telegraph Corp.) alcan- 
çar a posição do maior grupo de electrónica e telecomunicações do mundo. 


a Sfandard Elecfrica, SARL 
E UMA ASSOCIADA DA IT'T 
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O EMPREGO DO CIMENTO BRANCO... 


permite acabamentos 

mais perfeitos, mais 

duradouros e muito 
mais económicos 


Estude a vantagem do emprego do 


CIMENTO BRANCO LUSO 


IBRA 


Consulte os distribuidores gerais 


No Sul: SCIAL —T. do Corpo Santo, 15 — Telef. 66 61 86 — Lisboa 
No Norte: SCIAL — R. do Bonjardim, 205 — Telef. 25779 — Porto 


COLUNAS PARA 

ILUMINAÇÃO PÚBLICA 

MANILHAS PARA ESGOTOS 
dé 


SOCIEDADE 


PORTUGUESA 
E 


POSTES DE BETÃO PARA 
LINHAS ELÉCTRICAS 


MOSAICOS E AZULEJOS 


RUA D. ESTEFÂNIA, 94-A 
TELEF. 47812-50129 
LISBOA 
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COMPANHIA PORTUGUESA 
DE FORNOS ELÉCTRICOS 


PRODUTORA DE 


FERR QUO NS A 


FERRO LIGAS 


FERRO -MANGANES 
SÍLICO-MANGANÊS 


FERRO-SILÍCIOS 


15º) - 25% - 45% - 75º) - 90% 


* CARBONETO DE CÁLCIO 
CIANAMIDA CÁLCICA 
PASTA PARA ELECTRODOS 


E EM MONTAGEM 


GRAFITE 


FÁBRICAS SEDE 
CANAS DE SENHORIM Largo de S. Carlos, 4 2.º — LISBOA 
Telefone 67222 Telefones 325343-355 56-36 89 89 
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ISOLAMENTO 


o elulite TÉRMICO has 


CONSTRUÇÕES 


PAVIMENTOS 


E 


TECTOS 


* 


PAREDES 


* 


CAMARAS 
FRIGORÍFICAS 
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ISOLAMENTOS 
INDUSTRIAIS 


Terraços e pavimentos isolados com CELULITE no 
edifício da Comp.' de Seguros Fidelidade — P. das Águas Livres — Lisboa 


SETH, LDA. + RUA FILIPE FOLQUE, 10-1.º + TEL. 7301 56 — LISBOA 


| SAM DE PQ, 
Para e terraços, (EE esaLtico o 
paredes, fundações, etc. ip sa SS 


COLADEM DE 
TACOS DE MADEIRA 


TELEF. 730156  TELEG. EPALDA — LISBOA 
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Representantes em Portugal 
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EMPRESA TECNICA DE EQUIPAMENTOS ELÉCTRICOS, 


Rua Sá da Bandeira, 766-1.º E. 


Rua Rodrigo da Fonseca, 103, rc 


Tel. 24818 — PORTO 


LISBOA — Tel. 686072 
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Um perfil criado para 
a construçao moderna 


Edifícios administrativos, fábricas, estabelecimentos, hangares, 
pavilhões, hoteis e outras grandes construções erguem-se hoje 
segundo os princípios modernos da arquitectura, pelo processo 
de estruturas metálicas. Para vedar as vidraças, as placas dos 
parapeitos, etc., que se colocam nessas estruturas, obtem-se ex- 
celentes resultados com os perfis de vedação fabricados com a 
borracha de cloroprene Perbunan C. 

Estes perfis de vedação são muito flexíveis e vedam as vidraças 
ou outros elementos de construção, mesmo quando batidos por 
fortes terfpestades e chuvas. Em virtude das suas características 
de elasticidade, os perfis de Perbunan C, com uma composição 
adequada, possuem uma excelente resistência nas suas formas, 
quer dizer, mantém a sua forma primitiva durante longos anos, 
mesmo quando sujeitos a pressão e tracção. Resistem ao calor, 
de maneira que a sua função vedante não é prejudicada pelas 
temperaturas elevadas que se observam nas fachadas dos edifi- 
cios durante o verão. Possuem igualmente uma excelente resi- 
stência à acção da luz, do oxigénio, do ozono e dos fumos das 
fábricas. 
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E muito simples a montagem das vidraças e das chapas dos para- 
peitos com estes perfis. Aplica-se, tanto nos próprios perfis como 
nos caixilhos, uma cola de contacto fabricada igualmente à base 
de Perbunan C. Em seguida, aperta-se o perfil contra a peça, 
obtendo-se assim uma colagem segura. 

Os perfis de Perbunan C podem ainda utilizar-se na construção 
de barracões e de casa pré-fabricadas. Neste caso, as chapas 
ligam-se umas às outras por meio de perfis em | fabricados com 
Perbunan C (borracha de cloroprene). 

Escreva, por favor, aos nossos representantes, onde poderá obter 
informações mais detalhadas sobre os perfis de Perbunan C. 


S.A.R.L. QUIMICOR, Representante em Portugal 


Lisboa: Rua Sociedade Farmacêutica, 3, Tel. 42194 
Porto: Rua de Santos Pousada, 441, Tel. 54141 
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PORTO DA PRAIA DA VITÓRIA (AÇORES) 


ESTRAGOS NO QUEBRAMAR DEVIDO AO TEMPORAL DE 
26-27 DE DEZEMBRO DE 1962 


1— INTRODUÇÃO 


A baía da praia da Vitória está situada na 
costa este da ilha da Terceira, Arquipélago dos 
Açores. Este está localizado no Atlântico Norte, 
estando compreendido entre 36º 55' e 39º 43º de 
latitude norte e 24º 46" e 31º 16" de longitude 
oeste. A ilha Terceira faz parte do grupo central 
de ilhas do Arquipélago (figura 1). 

O esquema portuário que foi construído na 
baía da Praia da Vitória destina-se, fundamen- 
talmente, a permitir a descarga de combustíveis 
líquidos qualquer que seja o estado de agitação 
do mar à entrada da baía. 

O conjunto das obras marítimas é constituído 
por um quebramar de talude, enraizado na 
Ponta do Espírito Santo, e com cerca de 600 me- 
tros de comprimento, e por uma instalação de 
descarga, independente do molhe, composta por 
uma plataforma central e por duques de Alba. 
A plataforma central e os duques de Alba estão 
ligados entre sie a terra por meio de passadiços 
(figuras 1 e 3). 


2 — DESCRIÇÃO DA OBRA 


Dada a localização da baía da Praia da Vitória 
na costa este da ilha Terceira e a configura- 


por JOSÉ REIS DE CARVALHO 
Eng.º Civil (1.5. T.) 


Especialista do Laboratório N. de Engenharia Civil 


ção do contorno da ilha nessa zona, o quebramar 
só poderá ser atingido com violência apreciável 
pelos temporais de NE, E e SE. 

Note-se que para esta última direcção as ilhas 
de S. Miguel e Santa Maria constituem uma 
limitação à geração da agitação, reduzindo o 
«fetch» a um valor máximo de 120 milhas marí- 
timas. Este facto, aliado à configuração da bati- 
metria na extremidade sul da baía, faz com que 
os temporais de SE não possam atingir o que- 
bramar com grande violência (para um vento de 
40 nós a altura significativa da vaga em frente 
do quebramar é de apenas 4,00 metros). 

Por outro lado, a configuração batimétrica em 
frente da Ponta da Má Merenda, a norte da baía, 
determina uma forte dispersão de energia para 
os temporais de NE que, em face deste facto, 
atingem o quebramar com apenas cerca de 65 º/o 
da sua altura ao largo. 

Deste modo, só os temporais com direcções, 
em águas profundas, próximas de E podem atin- 
gir o quebramar com toda a sua violência. 

No estudo do dimensionamento do quebra- 
mar, a fixação da altura da onda de cálculo 
baseou-se nos elementos constantes da publica- 
ção «Sea and Swell Charts — North Atlantic 
Ocean», nos registos de velocidades e direcções 


* Comunicação apresentada à IX Conference on Coastal Engineering, Lisboa, 1964. 
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Fig. 1 — Porto da Praia da Vitória (Açores) 


de vento da estação meteorológica do aeroporto 
das Lajes (ilha Terceira) e nos valores das obser- 
vações locais sobre o estado de agitação do mar. 

Da análise comparada das duas primeiras fon- 
tes de informação chegou-se à conclusão que os 
ventos do quadrante de E ocorrem durante cerca 
de 20 º/o do tempo e não excedem a velocidade 
de 27 nós. 

Faz-se já aqui notar que, sendo os temporais de 
E gerados pelos ventos deste rumo mas soprando 
sobre a zona do Atlântico entre os Açores e o 
Continente Europeu, era nesta zona que interes- 
sava conhecer a velocidade máxima dos ventos 
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com aquela direcção e as dimensões da área 
sobre a qual ela poderia actuar. 

As observações locais indicavam que em No- 
vembro de 1955 tinham sido registadas vagas 
com 5,50 m durante um temporal de E e que em 
Fevereiro de 1956 foram estimadas vagas entre 
5,50 e 6,00 metros de altura. No entanto, o 
valor destas observações foi considerado, pelos 
projectistas, bastante reduzido, não tendo por 
isso servido para a fixação da onda de cálculo 
do quebramar. 

Tanto quanto se sabe, também não houve a 
preocupação de prever, a partir das cartas sinó- 


QuE BRAMAR 


ticas de tempo, quais teriam sido as alturas de 
vaga, que na realidade teriam atingido nesses 
períodos a baía da Praia da Vitória. 

A altura de vaga para o dimensionamento do 
quebramar foi determinada a partir do método 
de Sverdrup-Munk, utilizando os ábacos revistos 
por Bretschneider e considerando um vento com 
uma velocidade de 27 nós. A altura significativa 
é então de cerca de 5,20 metros, valor que foi 
fixado como o máximo para os temporais com 
direcções ao largo entre NE e SE. O período que 
se arbitrou para esta agitação foi de 12 segundos. 


3 — PERFIL DA OBRA 


O perfil dimensionado com base neste valor 
de altura máxima de vaga é formado por um 
núcleo de enrocamento de categoria C (até 
3,5 toneladas), revestido de enrocamentos da ca- 
tegoria B (3,5-5,5 toneladas) e da categoria A 
(12-13,5 toneladas), estes dispostos nas zonas do 
perfil directamente sujeitas à acção da vaga (fig. 2). 

Dada a dificuldade de se obter na pedreira a 
percentagem necessária de enrocamentos da ca- 
tegoria A, foram modificados os limites das ca- 
tegorias, permitindo-se que na camada interior da 
categoria A fosse utilizadó enrocamento de 9-12 
toneladas enquanto que na categoria B se permitiu 
o emprego de enrocamento de 2-3,5 toneladas na 
camada interior (figura 2). O peso específico dos 
enrocamentos é de 174,3 lbs'cuf (2,79 ton'm'). 

A aplicação da fórmula de Hudson ao perfil 
e para um valor da altura de agitação de 5,18 m, 
conduz a um coeficiente k de: 


RÃ so XE 
W(S.—1)' cota 12,7 x 1,72'x 2 


” 


3,00 


valor inferior ao aconselhado por Hudson o que 
leva a concluir que o binário, ângulo do talude- 
-peso dos enrocamentos, foi estabelecido com 
certa margem de segurança. 

Note-se contudo que o facto do enrocamento 
A ter sido dividido em duas categorias (A de 
12-13,55 e A — de 9-12 toneladas) diminuiu a 
segurança da obra, uma vez que a fuga de ele- 
mentos da camada exterior, sempre possível, se 
se verificasse, como veio a acontecer, a ocorrência 
de temporais com altura superior à de cálculo, 
expunha à acção directa da vaga elementos de 
peso inferior. Para esta camada o coeficiente k 
tem então o valor: 


O >< 12,7 
k = cds dai = 3,65 
10,43 


superior ao aconselhado por Hudson. 

Verifica-se portanto que, dada a alteração 
que foi introduzida nos enrocamentos da catego- 
ria A, a obra não tinha uma margem de segu- 
rança conveniente o que poderia ocasionar, no 
caso de um temporal mais violento que o de 
cálculo e de duração prolongada, um processo 
de destruição em cadeia com graves consequên- 
cias, uma vez que as camadas que sucessiva- 
mente ficavam expostas à fúria do mar têm uma 
estabilidade cada vez menor e portanto seriam 
destruídas, cada vez mais fâácilmente pelo tem- 
poral. Foi o que veio a acontecer nos fins de 
Dezembro de 1962. 

Pensando obviar os inconvenientes resultantes 
da subdivisão dos enrocamentos da categoria A, 
o projecto do molhe previu a colocação de blo- 
cos artificiais de forma especial, «cast through 
stones», englobando, pelo modo como eram colo- 
cados, as duas camadas do enrocamento. A colo- 
cação destes blocos artificiais, cujo peso é de 
14,5 toneladas, foi feita sem regras rígidas apenas 
obedecendo ao princípio de que pelo menos 1 
em cada 10 enrocamentos da categoria A, colocados 
no manto exterior deveria ser um «through stone». 

Note-se que, do ponto de vista da estabilidade 
do talude, o peso de cada «through stone» pouco 
ultrapassa o peso dos elementos da camada su- 
perficial, pelo que a sua capacidade individual 
de resistência pouco diferia dos referidos ele- 
mentos. Por outro lado o facto de cada bloco 
artificial estar rodeado de enrocamentos naturais 
faz com que a estabilidade do conjunto seja a 
dos elementos de menor capacidade de resistên- 
cia, tanto mais que a percentagem dos «through 
stone» é muito pequena. A sua presença, por- 
tanto, pouco ou nada contribuiu para a estabili- 
dade do talude. Foi aliás o que os ensaios em 
modelo reduzido, realizados para análise das 
avarias do molhe, mostraram. 

Três outros aspectos do perfil da obra, importa 
ainda focar. O primeiro diz respeito à cota 
inferior do manto exterior que no caso pre- 
sente se situa a um nível de 17'.0 abaixo do 
B. M. M. Como em baixamar de águas vivas o 
nível de água pode descer cerca de 1'.0 abaixo 
do B.M.M., a fixação do nível de (— 17'.0) não 
obedece estritamente ao indicado correntemente, 
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o que poderia levar a pensar que a ruína da obra 
teria começado exactamente nesta zona da pas- 
sagem dos enrocamentos da categoria A para a 
categoria B. (Os ensaios vieram mostrar que tal 
não se deve ter verificado. Aliás, de acordo com 
as experiências de Per Anders Hedar, as avarias 
num talude de enrocamento com uma inclinação 
de 2/1 não descem abaixo de um nível situado a 
0,8 H abaixo do B.M.A.V. o que está de acordo 
com os resultados dos ensaios realizados. 

Na face interior o manto de enrocamentos da 
categoria A estende-se até à cota ( + 5',0) ficando 
portanto a descoberto os enrocamentos da cate- 
goria B para todos os níveis de maré com excep- 
ção do praiamar de águas vivas, solução que não 
é correntemente indicada para obras em que se 
admite a existência de galgamentos. O «Technical 
Report n.º 4-—Shore Protection Planning and 
Design» aconselha que os enrocamentos do manto 
exterior desçam, na face interior, até ao 
B.M.A.V. 

Mas foi ma cabeça do molhe que o projecto 
mais se afastou das regras e indicações normal- 
mente aceites. Com efeito é do conhecimento 
geral que as solicitações que se verificam nos 
elementos que constituem a cabeça de um molhe 
do tipo do quebramar da Praia da Vitória, são 
superiores às que se verificam nos elementos do 
perfil corrente dado o jacto que a rebentação da 
vaga ocasiona nesta zona. Por outro lado a cur- 
vatura do molhe, na zona da cabeça, determina 
que a estabilidade dos elementos que constituem 
o quadrante interior da cabeça, seja inferior à 
dos elementos do perfil corrente. A conjugação 
destes dois factores faz com que a cabeça dos 
molhes seja uma zona em que se torna necessário 
utilizar elementos de maior peso sobretudo no 
sector interior e na zona compreendida entre o 
mínimo e o máximo nível de água. O «Technical 
Report n.º 4» indica, para o caso de taludes de 
enrocamento, um acréscimo de peso da ordem de 
10 º/o na zona da cabeça. 

Daqui se conclui que a cabeça do molhe da 
Praia da Vitória não tinha o mesmo coeficiente 
de segurança que o perfil corrente. E no entanto, 
e por razões que se indicarão adiante, a cabeça 
do molhe apenas teve ligeiras avarias no sector 
interior. Faz-se contudo notar que a diminuição 
de resistência da cabeça de um molhe é alta- 
mente perigosa uma vez que qualquer avaria que 
aí se verifique se estende facilmente ao perfil 


corrente mesmo que este, por si só, possa resistir 
ao temporal. 


4— O TEMPORAL DE 26-27 DE DEZEMBRO 
DE 1962 


Nos últimos dias do mês de Dezembro de 
1962 foi o arquipélago dos Açores assolado por 
forte temporal. Os elementos meteorológicos 
cedidos pelo Serviço Meteorológico Nacional 
(figura 4) mostram que a partir das 00,00 h do 
dia 25 de Dezembro uma depressão alongada, 
centrada sobre a ilha de Santa Maria, ocasionava 
no seu bordo norte, fortes ventos de E. Esta 
depressão manteve-se sensivelmente estacionária, 
por vezes com mais de um núcleo, até depois do 
dia 27, embora a partir desta data tenha dimi- 
nuído de actividade. 

No período de maior intensidade o vale depres- 
sionário estendia-se dos Açores ao Norte de 
África originando ventos de E de 25 a 40 nós 
desde a costa de Portugal Continental aos Açores. 

Estas condições, particularmente desfavoráveis, 
determinaram uma área de geração suficiente- 
mente extensa e uma duração de vento que per- 
mitia, para velocidades de vento até 40 nós, a 
geração do estado de completo crescimento da 
agitação. 

Infelizmente os dados meteorológicos existentes 
não permitiram definir perfeitamente o campo 
das velocidades na zona de geração. Apenas foi 
possível limitar a velocidade do vento que oca- 
sionou o temporal, entre 30 e 34 nós, o que se 
traduz por uma altura significativa da agitação, 
que pode ter variado entre 6,50 metros (30 nós) e 
9,00 metros (34 nós) de acordo com os gráficos 
de Neumann. Faz-se notar contudo que a veloci- 
dade do vento de E na estação meteorológica das 
Lages (Ilha Terceira) não ultrapassou 25 nós. 

Segundo os elementos colhidos junto da Junta 
Autónoma do Porto de Angra do Heroísmo (Ilha 
da Terceira), o quebramar da Praia da Vitória foi 
continuamente galgado desde a tarde do dia 25 
de Dezembro até à manhã do dia 28, com espe- 
cial relevo nos períodos de preiamar. Nestes 
períodos a «gerbe» que resultava do galgamento 
do molhe passava sobre a própria estrutura de 
atracação dos petroleiros. 

Durante os dias 26 e 27 foram estimadas as 
alturas da agitação por meio de um método espe- 
dito. Os valores encontrados, média dos registos 
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feitos durante períodos de 30 a 40 minutos, 
foram os seguintes: 


Data | ie H Data Rm ) H 
a E E | 
26-12-62| 1000 |- 8.c0m ||27-12-62]| 10.00 |= 8.7om 

, 14.00 | = 8.10M » 15.00 | = 8.00m 
» | 18.00 |= 8.7om » 15.00 |= 87om 


Os registos mostraram que as variações das 
alturas individuais não ultrapassaram, em cada 
período de observação, estes valores médios em 
mais de 0,90 m. 

Conjugando estes valores com a análise, tão 
pormenorizada quanto possível, dos campos de 
velocidades que se devem ter verificado sobre o 
Atlântico Norte chega-se à conclusão que o tem- 
poral que assolou o porto da Praia da Vitória 
deve ter tido início às 00,00 h (T.M. G.) do dia 
25 de Dezembro. Às 12,00h (T.M.G.) já a 
altura significativa da agitação tinha atingido 
cerca de 3,00 metros para às 00,00 h (T.M. G.) 
do dia 26 de Dezembro atingir 6,00 metros. 
O temporal, no entanto, ainda se tornou mais 
violento durante este dia, devendo ter atingido a 
sua máxima violência cerca das 18,00 h (T.M.G.). 
A partir deste instante o estado de agitação do 
mar manteve-se sensivelmente constante até 
cerca das 18,00 h (T.M.G.) do dia 27, dimi- 
nuindo depois rápidamente de intensidade. O 
período associado ao temporal foi, nos períodos 
de maior violência, de 13 a 14 segundos. 


5 — ESTRAGOS NO QUEBRAMAR, 
SUA ANÁLISE 


Na altura em que se verificou o temporal o 
quebramar estava quase totalmente construído, 
faltando apenas colocar os enrocamentos da ca- 
tegoria A entre as cotas (+ 12'.0) e (+-22'.0) 
entre o enraizamento (perfil 0) e o perfil 175 m. 

Os primeiros sinais de avarias foram obser- 
vados na manhã do dia 26 de Dezembro quando 
se verificou, visualmente, que algumas pedras 
tinham desaparecido na zona emersa do perfil 
335 metros, onde os enrocamentos da categoria A 
continuavam a ser desalojados das suas posições 
e rolavam ao longo do perfil desaparecendo 
depois debaixo de água, sendo mais tarde encon- 
trados, nos levantamentos efectuados após a 
tempestade, no pé do talude. 
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Pode-se daqui já extrair uma primeira conclu- 
são. A ruína do molhe foi devida à falta de 
estabilidade do talude exterior e não devido aos 
galgamentos anormais do molhe. Estes teriam 
determinado um tipo de ruína que começaria 
pelo deslocamento dos blocos do talude interior. 

Na noite de 26 para 27 de Dezembro, tendo 
o temporal atingido a sua máxima violência, foi 
o molhe particularmente atingido e avariado. As 
avarias generalizaram-se e na manhã do dia 27 
já era visível uma profunda brecha imediata- 
mente junto à cabeça, onde o temporal cortou 
todo o perfil acima da cota (+ 0.00) erodindo 
inclusivamente os enrocamentos da categoria € 
(núcleo). 

Durante o dia 27 o molhe continuou a ser 
destruído a ponto dos galgamentos, que então se 
verificavam sobretudo no preiamar, terem posto 
em risco a própria instalação de atracação, aliás 
avariada já durante a noite anterior. Felizmente 
durante a noite de 27 para 28 o temporal dimi- 
nuíu de violência de modo que no dia 28 não 
se observaram quaisquer novas destruições. 

O balanço do temporal pode resumir-se da 
maneira seguinte (figuras 3, 5,6 e 7): 


— perfis O a 176 metros: neste troço do mo- 
lhe, completo até à cota (+ 12'.0) apenas se 
verificaram ligeiras avarias superficiais; por 
outro lado o talude interior e exterior ficou 
revestido por um volume considerável de 
pedras de pequenas dimensões que foram 
transportadas da costa a norte do molhe; 

— perfis 176 a 291 metros: este troço do mo- 
lhe ficou em boas condições tendo apenas 
fugido alguns enrocamentos da categoria A 
do talude exterior e assentado ligeiramente 
o seu coroamento; contudo um número 
apreciável de enrocamentos da categoria B 
foram deslocados no talude interior, abaixo 
da cota (+ 5'.00), o que leva à conclusão de 
que os galgamentos tiveram neste troço 
efeito bastante pernicioso ; 

— perfil 291 a 442 metros: este troço foi um 
dos mais duramente atingidos pelo tempo- 
ral, tendo ruído todos os enrocamentos da 
categoria A, tanto no talude exterior como 
no coroamento, elementos que rolaram pelo 
talude exterior, juntamente com pedras da 
categoria B que estavam abaixo dos enroca- 
mentos da categoria A, para o pé do talude; 


Fig. 3 — Vista geral do porto 


no entanto conservaram-se ainda em posi- 
ção, no talude interior, alguns elementos da 
categoria A e «through stones» se bem que 
em estabilidade muito precária; se o tem- 
poral tem tido uma duração um pouco maior 
estes elementos teriam também ruído e este 
troço do molhe teria então ficado arrazado 
a uma cota vizinha de (+ 0.00); o tipo 
de avaria verificada nesta zona do quebra- 
mar confirma o que se verificou nos troços 
anteriores e indica que o processo de ruína 
se iniciou pelo talude exterior ; 

— perfis 442 a 530 metros: este troço foi menos 
atingido que o anterior pois aproximada- 
mente 45 metros ficaram quase intactos; 

— perfis 530 a 565 metros: foi este o troço 
em que se registaram as maiores destrui- 
ções; o quebramar foi arrazado até à cota 
(+ 0.00); 

— perfis 565 a 585 metros (cabeça do molhe) : 


este troço ficou em boas condições apenas 
se tendo verificado a fuga de alguns enro- 
camentos, na zona interior e abaixo da li- 
nha de água, confirmando assim o conhe- 
cido sobre o processo de ruína das cabeças 
dos molhes. 


A primeira conclusão que se extrai das ava- 
rias verificadas ao longo do molhe é que o 
processo de ruína foi extremamente irregular; 
a cabeça ficou quase intacta, brecha no troço 
consecutivo, depois cerca de 90 metros com pe- 
quenas avarias para logo a seguir, e numa exten- 
são de 150 metros, o molhe ter sido duramente 
avariado. 

Esta extrema irregularidade teve contudo 
a vantagem de permitir deduzir, na parte 
visível, o processo de ruína do molhe, bem evi- 
dente nos vários graus de ruína que se indica- 
ram anteriormente: inicialmente foram os ele- 


TRONICA 
569 


“g 

a R 

[sm 

NE 

a 

s 

E” = 8 
: < 
ç n 
AN £ 
n 

q 

: o 

s 


EISLICIAR 


MONTAGEM EM BRANCO DUMA BOMBA DE ACUMULAÇÃO HPV 270, 
17 m3/seg., 40.000 CV, DESTINADA À CENTRAL DO ALTO RABAGÃO 


E CONSTRUÍDA POR CONSTRUÇÕES METALOMECÂNICAS MAGUE, 
S.A.R.L., SEGUNDO PROJECTO E TÉCNICA SULZER FRÊRES. 


PONTES ROLANTES ———— —— — Projecto e fabrico 

GUINDASTES — — Projecto e fabrico 

APAR. DE ELEVAÇÃO ESPECIAIS — Projecto e fabrico 
TURBINAS HIDRÁULICAS —— Fabrico segundo licença de A.C. M. de Vevey 
TURBINAS A VAPOR ——— > Fabrico segundo licença de Brown Boveri, Cie 


EQUIPAMENTOS INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGQUE ;.-. 


ALVERCA DO RIBATEJO- PORTUGAL 


TECNICA XIII 


aços vazados 


ferros meelhanite 
mecânica mévia E PESADA 
construções metálicas 


COMPORTAS 
TURBINAS , 


CONDUTAS FORÇADAS 
RESERVATÓRIOS SOB PRESSÃO 
PERMUTADORES DE CALOR 


CaRSa +15 À . 
” ad go qi 4 
X E OP e ' pe 
we 


Centrais Hidráulicas 


Válvula Esférica para Lomaum 


DIVISÃO METALO-MECÂNICA 
COMPANHIA UNIÃO FABRIL 


AVENIDA 24 DE JULHO, 170 — LISBOA 2 4 TELEFONES 670421 E 672061 + TELEGRAMAS: METALOFABRIL 


TECNICA XIV 


Fig. 5 — Vista do temporal e das avarias sofridas pelo molhe 


mentos da categoria A que foram arrastados ro- 
lando pelo talude abaixo; depois foram os enro- 
camentos da categoria B que, já sem protecção, 
foram deslocados rolando também para o pé do 
talude; finalmente as vagas após terem também 
arrastado elementos do núcleo, empurraram os 
poucos elementos do talude interior, ainda em 
posição, para o lado do porto abrindo uma brecha 
no molhe como se verificou junto à cabeça do 


molhe. Não foi contudo possível a partir da 
análise dos levantamentos efectuados em Janeiro 
de 1963 deduzir o processo de ruína nas zonas 
submersas do perfil uma vez que aí se acumu- 
laram os enrocamentos que foram deslocados 
na parte superior do talude. 

Na parte interior do perfil os galgamentos 
apenas deslocaram alguns enrocamentos da cate- 
goria B. 
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Dois factores poderiam ter determinado esta 
irregularidade do processo de ruína: variações 
pronunciadas da altura da vaga ao longo do 
quebramar ou então irregularidades na cons- 
trução que o temporal, superior ao de cálculo, 
teria posto em evidência. 

O traçado de dois planos de ondulação, um 
com a direcção E ao largo e outro com a direc- 
ção E— 10º — N, mostrou que, para o primeiro 
rumo, o ataque da onda se fazia frontalmente 
verificando-se as seguintes variações de altura ao 
longo do quebramar: ligeira concentração no 
enraizamento, depois ligeira diminuição para a 
zona da cabeça onde, localmente, se verificou um 
abaixamento sensível, embora no troço imediata- 
mente antes houvesse uma ligeira concentração 
de energia. 

Para o segundo rumo a vaga atacava a obra 
com um ângulo pequeno, não indicando qualquer 
variação da altura da vaga ao longo do quebra- 
mar. A análise das fotografias feitas no decorrer 
do temporal mostra que o ataque se fez sempre 
frontalmente. 

A variação da altura da vaga ao longo do 
molhe para a vaga de E ao largo explica a não 
existência de avarias na testa e a brecha no 
perfil adjacente, mas não consegue explicar o 
estado em que ficou a parte restante da obra 
onde alternam troços com avarias graves com 
troços praticamente intactos. Independentemente 
do facto do quebramar ter sido atingido por 
vagas com altura superior à de cálculo e durante 
um largo espaço de tempo, aquelas diferenças 
parecem só serem explicadas por diferenças cons- 
trutivas. 


6 — ENSAIOS EM MODELO REDUZIDO 


Com o intuito de esclarecer alguns pontos do 
processo de ruína do quebramar do porto a 
Praia da Vitória considerou-se de muito interesse 
realizar alguns ensaios em modelo reduzido. 

Estes decorreram em canal de vagas conside- 
rando-se portanto apenas o caso do perfil cor- 
rente atacado frontalmente pela vaga. 

Os ensaios, realizados à escala 1/50, foram 
inicialmente conduzidos para os níveis do preiamar 
e baixamar, sem variação de nível, e finalmente 
com variação de nível de modo a reproduzir a 
maré local, 

Mostraram os ensaios que a reprodução da 
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maré, ainda que por saltos, conduz a maiores 
avarias, para a mesma altura de vaga, que os 
realizados com um único nível de água, Por este 
facto as experiências foram conduzidas na sua 
parte final, com reprodução da maré e para várias 
alturas máximas de vaga. 

Todos estes ensaios mostraram que as primei- 
ras avarias ocorreram para vagas de cerca de 
6 metros o que mostra que o molhe estava bem 
concebido tendo em atenção a altura de vaga 
utilizada no cálculo. 

AS primeiras avarias ocorreram na zona entre 
os níveis de baixamar e preiamar, rolando as 
pedras ao longo do talude e caindo na sua base. 
Depois, à medida que a altura da vaga ia cres- 
cendo a zona em que se registavam fugas do enro- 
camento era cada vez maior, estendendo-se tanto 
para o coroamento do molhe como para a zona 
submersa. Note-se contudo que os enrocamentos 
da categoria B abaixo da cota (— 17'.0) nunca 
foram deslocados, se bem que para as alturas de 
vaga mais elevadas esta zona já estivesse reves- 
tida com enrocamentos da categoria A que tinham 
sido deslocados da parte superior do talude. 

O coroamento do molhe foi atingido, no pro- 
cesso de ruína, para altura de vaga de 7,50 metros. 

Conclui-se daqui que a altura máxima do tem- 
poral foi sensivelmente igual a este valor, pelo 
que os ensaios prosseguiram conservando esta 
altura durante o tempo correspondente aquele 
que em Praia da Vitória se mantiveram as con- 
dições mais desfavoráveis da vaga. 

Note-se que no paramento interior do molhe, 
os enrocamentos da categoria B abaixo da cota 
(+-5'.00) começaram a ser deslocados para altu- 
ras de agitação de cerca de 7,00 metros, o 
que vem corroborar a afirmação de que a altura 
máxima do temporal foi de cerca de 7,50 metros. 
Com efeito, para alturas superiores as avarias 
que se verificaram no talude interior seriam 
muito superiores do que as realmente se obser- 
varam em Praia da Vitória. 

Também no talude exterior um aumento 
mesmo pequeno na altura de vaga acelerava 
extraordinâriamente o processo de ruína deter- 
minando, para alturas de vaga de 7,70 m, a forma- 
ção de uma brecha do tipo da que se abriu junto 
à cabeça do molhe. Esta brecha pode, portanto, 
ser perfeitamente explicada pela ligeira concen- 
tração de energia que os planos de ondulação 
mostraram verificar-se nessa zona do molhe. 


Se bem que para a onda máxima de 7,50 metros 
as avarias verificadas nos vários ensaios realizados 
fossem ligeiramente diferentes, não foi, con- 
tudo, possível encontrar uma explicação para a 
grande irregularidade das avarias que só se 
baseia em possíveis irregularidades na constru- 
ção do molhe, possivelmente resultantes de, em 
certos troços, terem sido colocados enrocamen- 
tos da categoria A de peso superior ao previsto 
no projecto. Com efeito, segundo informações 
colhidas localmente teve-se conhecimento da 
colocação de enrocamentos de 20 toneladas, em- 
bora não tenha sido possível saber em que 
perfis tais enrocamentos foram colocados. 


7 — CONCLUSÕES 


A análise conjunta da observação das avarias 
verificadas no quebramar da Praia da Vitória e 
dos ensaios em modelo reduzido permite extrair 
as seguintes conclusões : 


-— deve haver particular cuidado na determina- 
nação dos maiores temporais que podem 
assolar uma obra marítima; uma deficiente 
fixação desse valor pode ocasionar sérias 
avarias como no caso presente; deve por 
isso haver a preocupação de estabelecer 
programas de observação com equipamento 
conveniente; além da determinação da pro- 
babilidade da ocorrência de vagas com deter- 
minada altura, importa ainda conhecer o 
máximo tempo possível de duração conse- 
cutiva desses temporais (no caso da Praia 
da Vitória a duração prolongada do tempo- 


ral foi um dos factores determinantes do 
grau de avaria); 

— todo o projecto de quebramar deve ser su- 
jeito a ensaios em modelo reduzido não só 
para averiguar qualquer deficiência de pro- 
jecto mas também conhecer o grau de segu- 
rança da obra e o seu comportamento para 
alturas de onda superiores à de cálculo; 

— tendo-se verificado um paralelismo muito 
aceitável entre as avarias observadas na na- 
tureza e nos ensaios em modelo, estes per- 
mitem, em conjugação com o conhecimento 
do regime de agitação junto a um quebramar, 
o seu dimensionamento de modó que o 
custo total seja mínimo. 


O autor exprime os seus agradecimentos ao 
engenheiro Leiria Gomes, director da Junta Autó- 
noma do Porto de Angra do Heroísmo, que tor- 
nou possível, com o envio dos elementos de 
campo necessários, a realização deste trabalho. 
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SUMÁRIO 


O quebramar do porto da Praia da Vitória (Açores), sofreu graves avarias durante 


o temporal de 26-27 de Dezembro de 1962. 


O autor analisa as avarias e compara-as com as obtidas em ensaio em modelo 
reduzido. Conclui sobre a premente necessidade de se estabelecer um programa sério 
de observações do regime de agitação nas zonas onde se pretende construir portos 
ou ampliar os já existentes e sobre a utilidade dos ensaios de estabilidade, das obras 
marítimas, apontando a concordância perfeita entre as avarias sofridas pelo quebramar 
da Praia da Vitória e as observadas em modelo. 
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SUMMARY 


The breakwater of Praia da Vitória harbour (Azores) was seriously damaged 
during the storm of 25th-2ath December 1962. 

The author analyses the damages, comparing them with the results of model 
tests, and concludes that comprehensive programs of observations of wave conditions 
are urgently needed in areas where the construction or extension of ports is envisaged. 
The advantages of stability tests of maritime structures are indicated, the remarkable 
agreement between the damages underwent by the breakwater in Praia da Vitória and 
the results of the model tests being underlined. 


CONGRESSO SOBRE FUNDAÇÕES PROFUNDAS 


A Sociedade Mexicana de Mecânica do Solo com a colaboração da Sociedade Americana 
de Engenheiros Civis (ASCE), Divisão de Mecânica do Solo e Fundações, e a Sociedade Americana 
para Ensaios e Materiais dos Estados Unidos (ASTM), vai realizar de 7 a 12 de Dezembro de 1964 
um Congresso sobre fundações profundas, cujas línguas oficiais serão o espanhol e o inglês. 

Esta reunião, que terá lugar na cidade do México versará os seguintes assuntos : 


| — Análise teórica da capacidade de carga e dos movimentos horizontais e verticais de funda- 
ções profundas, comparada com valores de ensaios de campo. 

IH — Ensaios em modelos relacionados com o tema geral do congresso. 

HI — Ensaios e observações de campo relativas a fundações profundas. 

IV — Métodos e equipamentos usados na execução de estacas. 

V — Casos especiais de estruturas com fundação profunda. 


As comunicações deverão ser apresentadas até 15-10-64. 
Qualquer pedido de informações deverá ser digido a Prof. Rafael Arias Nuvillo — Comisión 
de Congresos y Reuniones — Vértiz 800, 5.º Piso — México 12, D. F. — México. 


TRONICA 
576 


PoclainT.Y.45 


Escavadora=-Rctro Escavadora-Grua 


e POTÊNCIA AUMENTADA 

e ROTAÇÃO TOTAL E CONTÍNUA 
e ESTABILIDADE 100% 

e 20 BALDES DIFERENTES 


mm 


REPRESENTANTES EXCLUSIVOS 

Sociedade de Mecanização Industrial e Agricola, S.A R.L. 
Avenida Padre Manuel da Nóbrega, 8-B.1.º - LISBOA 
Telefone: 724053/4/5 


TECNICA XV 


FANAF EL 


FÁBRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 


APART.: 9 OVAR ESCRITÓRIO E FÁBRICA 
TELEF.: 52093-52094 ESTRADA DE S. JOÃO 


TELEG.: FELTROS PORTUGAL 


FELTROS, FILTROS E TECIDOS INDUSTRIAIS PARA TODOS OS FINS: 


EM PEÇAS, PANOS, MANGAS OU SACOS, 
DE QUALQUER FIBRA 
PASTA PARA PAPEL TEGCELAGENS 


AS DE Ie SR “ACABAMENTOS 
FIBROCGIMENTO TINTURARIAS 


CURTUMES LAVANDARIAS 
FIAÇÕES DE ALGODÃO ESTAMPARIAS 
Wes Pai Wi E DS ESSE DA 
INDÚSTRIAS QUÍMICAS DESPOEIRAGEM 


IMPERMEABILIZAÇÕES 


TORNA TUDO ESTANQUE 


Impermeabilização e vedação perfeita de: 


Terraços, Empenas, Algerozes, Clarabóias, 
Paredes húmidas, Coberturas, Fendas, 
Juntas de dilatação, etc., etc. 


Fácil aplicação — Elasticidade excepcional 
Não endurece nem apodrece 


À base de borracha + alumínio 


Reforça a acção impermeabilizante do «Rubio», dando às construções um aspecto agra- 
dável — Protege todos os materiais contra a humidade, a ferrugem e a corrosão 


FLUIDO E MASTIC 


RUBIO E ROOFATEX dois produtos que se completam 


Representante exclusivo : 4 Cnlundo Femmes 


Rus Palmira, 46-A-B-C — LISBOA-1 — Telefs. 8301 63-84 26 39 — End. toleg. «ORFER» 


TECNICA XVI 


C. D. U. 534.144 


VIBRAÇÕES FORÇADAS DE SISTEMAS NÃO LINEARES 


COM 


1 — INTRODUÇÃO 


Num trabalho anterior [1] (*) abordámos o 
estudo das vibrações de sistemas não amorte- 
cidos de elasticidade do tipo k (x+texº), em 
regime livre, e com um grau de liberdade. 
No intuito de completar aquele trabalho propo- 
mo-nos agora apresentar alguns métodos para o 
estudo do referido sistema em regime forçado, 
tendo o cuidado de descrever sumariamente 
aqueles que são vulgarmente conhecidos para 
podermos efectuar uma crítica comparativa. 


2 — SISTEMAS COM ELASTICIDADE 
SIMÉTRICA 


O modelo matemático que consideraremos terá 
a forma 


mxtHk(x+tHex)=F(t) czo 


ou 
x + og (x + ex) =p (t) (1) 


e dº x 
em que x significa como habitualmente qe e 


sendo: 


m — massa pontual equivalente ; 

k — característica da mola linear equiva- 
lente nas zonas das deformações infi- 
nitesimais ; 

: — factor de não linearidade; 

)y — frequência própria angular (pulsação) 
do sistema linearizado em torno da 
posição de repouso; 


(*) Números entre parêntesis rectos referem-se à biblio- 
grafia indicada. 


ELASTICIDADE DO TIPO k(x+:x) 


Pelo Eng.* Mec. (1. S. T.) 


JORGE NEVES DA SILVA 


(Assistente no 1. 5. T.) 


Xo — amplitude da deformação x sempre 
positiva (xo > 0); 

F(t),F(vt) — função perturbadora periódica de fre- 
quência angular fundamental »v, de 
amplitude Fo; 

p(t),p(t) — função perturbadora normalizada 
(m==1) de amplitude p,. 


Introduzindo os coeficientes adimensionais 


r=Vex ; to=Vixo para :>0 
r=V=:* ; Xo == V-E xo para :< O 
Ve Ve Po 
Â — E ; Ao = 9 para E >: 0) 
tdo o 
Vs V =| 
á P, Ao = — e para :<O 
Do 1) 
v f, 
e pm mens 
Am fo 
J 
em que =— e f;= ne 
2”. 2% 
(1) ficará o 
À trt$=A (2) 
OM 


+, 


vindo os sinais + ou — consoante é : >=0 ou 


e O. 


Dado que o princípio da sobreposição não é 
válido para os sistemas não lineares, não tem 
interesse especial a consideração de uma pertur- 
bação sinusoidal pura. No entanto, porque, de 
certo modo, é uma perturbação padrão, normal- 
mente suporemos os sistemas actuados por fun- 
ções posinvt ou po cosvt. 


2.1. — Função perturbadora elíptica 


Este método tem a vantagem de ser exacto 
mas implica uma forma especial para a pertur- 
bação, sinusoide ou cosinusoide elípticas de 
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parâmetro variável com a amplificação, com o 
factor £ e com a relação de frequências. Hsu [2] 
designa esta perturbação ideal por «favorável», 
na medida em que é a função que se ajusta ao 
movimento, isto é, supõe-se que a função per- 
turbadora é tal que tem sempre a forma exacta 
do movimento do sistema perturbado, oscilando 
este em fase ou em oposição mas nunca desfa- 
sado (por não haver amortecimento). Transfor- 
ma-se assim um sistema forçado num outro auto- 
-excitado. 
Seja então na equação (2) 


Ap)=a x (3) 


em que « é uma constante a determinar. Em par- 
ticular 


Ao == à X0 
e (2) virá 
x A 
E +(1-D) p= (4) 
1) N Xo / 


A integração desta equação diferencial origina 
soluções dependentes do sinal do coeficiente 


A ES 
(1 — —e ) e do sinal de xº e pode efectuar-se 
xo , 


tal como fizemos em [1], embora, para abreviar 
a exposição, nos limitemos a verificar soluções 
propostas. 


virá positivo ou negativo consoante o sis- 
to 

tema vibra em fase ou em oposição de fase com a 

força perturbadora. 

Convirá apresentar os diagramas da deforma- 
ção elástica imaginária, e de fase, do sistema li- 
vre que caracteriza o modelo matemático obtido 
em (4). 

Uma vez que o sistema deve oscilar com uma 
frequência imposta, entre as quatro hipóteses 
admissíveis (fig. 1 a 4) só é possível encontrar 
soluções para à, << xo com : > 0 (fig. 1) ou 
E< O (fig. 2), e para À > xo com: >o0 (fig. 3). 
Nesta última hipótese o sistema oscilará em 
torno de um dos centros de estabilidade (a,o) 
ou (— a,o) desde que 


xo < | | e 
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Fig. 2 


ou em torno do ponto de equilíbrio para o sis- 
tema em repouso (0,0), desde que 


E >0 
E X e 
o> | x| Pt 


em que x; representa um dos pontos de fronteira 
(— x, 0) ou (+ x1,0) (ver fig. 3). 


Fig. 3 


Fig. 4 


Os centros de estabilidade obtém-se fazendo 


Fe=F()=(1— fo Jata=o 


Xo 


x as 
u=a=/ | 
Xo 


donde 


= gem / do — 1 
Xo 


Os pontos de tronteira virão dados por 
fé Fdr=o 


e, substituindo, 
Y + 2%w—-246=0 
ou 


Ão == 1 fãs x 


2.1.1. Solução em cn 


Suponhamos que uma solução da equação 
diferencial (4) é o coseno elíptico 


maos Cult 


em que o argumento u se relaciona com a fre- 
quência angular fundamental da perturbação (») 


' E 2K 
através da expressão u=-"— v (t — to) obtida 


T 
por desenvolvimento em série de Fourier da 
função cn u [3, pág. 13](*. K é o i.e.j. completo 
de 1.2 espécie. 


Derivando (5) duas vezes 


 ,sn(ulk)dn (u|p) 


tt == to 


e [2K 
ig | 
- 


+ sn? (u|2) en (u/p)] 


) [-cn (u/g) dn? (u |) + 


a... > E cos | (2s+1)zu 

N 1/2 pr Fm qês+ 2K 

com 
Tr 
E = fa E ae 
e Jo  |1—usin?0]l2 
x 
K'(u) = K(t —u) q=e. “K 
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e tendo em atenção as relações entre sn, cn e 
dn [3, pág. 16] 


, 3 (6) 
re to( 2 ) (29-1) en (ulp)-2p crê (u|8)] 


T 


Substituindo (6) e (5) em (4) e tendo em 
conta os coeficientes adimensionais, virá 


to) 241) (25 ) E (1 2) fentuja) + 


+ xo | 2u Fa ) ch + | cnº (u | H) ==() (7) 


wo) 


Para que a equação (7) se verifique indepen- 
dentemente dos valores da função cn (u's) é ne- 
cessário e suficiente que 


Te 
| (8) 
so Rd C+r=o0 
Te 
que pode tomar a forma 
( +), 
“a ES 
| 21-28 + 2x) 
K == K (1) Ven ntaça [3, pág. 106 109]) (8”) 
+ xó 


É eqato 


Cc = = 
2u (28) 
| T 


O parâmetro | só pode tomar valores entre 
zero (cn cos) e 1 (cn ch). Para que (5) seja 
uma solução periódica terá que ser c? — 0 e só 
serve o sinal + nas expressões (8'), e portanto 
dessas mesmas expressões se conclui que 


E 20 
A, << É guia 
o %o + Xo 
2 
O sistema oscilará: 
— para À, < xo, como se indica na fig. 1; 


| 1 
— para wi Lx + E xo”, como se in- 
dica na fig. 3, em torno da origem. 
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2.1.2. Solução em sn 


Suponhamos agora para solução de (4) o seno 
elíptico 
x==ro sn (uÊ) (9) 


em que o argumento u se relaciona com a fre- 
quência angular fundamental da perturbação (») 
pela expressão já apresentada em 2.1.1 


2 v (t—to) 


= 


do mesmo modo obtida através do desenvolvi- 
mento em série de Fourier de sn u [3, pág. 13] **º. 
Derivando (9) duas vezes virá como em (6). 


x = to(2S ) [1-5 sn (u/2) + 


1. 


+ 2: snº (u/9 | (10) 


e (7) será agora 


“| (2-2) ali " (1 - =) [sn (u/2) + 


+ m[24(20) ex: | snº (u/5) =0 (11) 


e tal como em (8) 


| (=E —31) E) + (1 — = j=o 


-(/2KN2 ,, (12) 
| 2: ( E) d+ti=o0 
podendo tomar a forma 
» RA % 
“O za + Mo 
K =2K (:) (tabelas, [3, pág. 106/109]) (12) 
c irao “E T; 
é o 
(o HG) 
' o 2- asi 4 + 1a (25+1) = 
() sn mad dk SEI E uma 
te K s=o | — q 2 K 
com 


dis 


K()=K (1—&8) q=e K 


O parâmetro É só pode tomar valores entre 
zero (sn =sin) e 1 (sn = sh). Portanto para 
que (9) seja uma solução periódica terá que ser 


cC>0 e só serve o sinal — das expressões 
(12). De (12') conclui-se que 
E <oOo 


do po Si (1 sos 2) 
O sistema oscila como se indicou na fig. 2. 


2.13. Solução em dn 
Supondo agora para solução de (4) a função 


x =x, dn (u/n) (13) 
em que u se relaciona com v através da ex- 
pressão 

K 


SA ” (t — to) 


obtida pela apreciação do desenvolvimento em 
série de Fourier de dn u [3, pág. 13] '*). 
Derivando (13) duas vezes 


ba (e ) [2-0 dn (uid) — 2 do? (u/o) | 


(14) 
e, tal como em (7), será agora 


Xo [en(S'e + (1 rd )| dn (u |n) + 
13 xy 


Ê K Y* 
+5|—2(5) ctg Jd (u |n) =0 
Tr - 


vindo em vez de (8) 


K '2 (16) 
| -2( + C+j= 0 
E 27 O q” sru 
(*) dn (u/n) = 4 2 sz cos ——. 
2K K qui 1-0" 
com 
, 
koto [mat 
"O [1 —n sin?6] /2 
mem: DS 
K()=K(1=-n) q=e K 


podendo escrever-se ainda 


2x +2 (xy — 45) 
em SE pe 
x 


K =K (n) (16”) 


O parâmetro n só pode tomar valores entre 
zero (dn==1) e 1 (dn ==sech), e então para que 
(13) seja uma solução periódica terá que ser 
cC>0 e só serve o sinal + em (16). Destas 
expressões se conclui que 


20 


oscilando o sistema como se indicou na fig. 3 
em torno dos pontos (— a,0) ou (+ a,o). 


2.1.4. Representação gráfica 


O gráfico I apresenta as famílias de curvas de 
amplificação do movimento do sistema, paramé- 
tricas em Ao — amplitude da força perturbadora. 

Os sistemas macios (=: 0) estão representados 
em ordenadas negativas e o movimento é do 
tipo definido pela equação (9). 

Os sistemas duros (: 0) estão representados 
em ordenadas positivas e o movimento é do tipo 
definido pela equação (5) no domínio 4, < x, + 


i a ê 1 s 
A E x). No domínio x; + > o My < xx 


o movimento existe também com a forma (13) 
desde que, evidentemente, seja esta a forma da 
perturbação (curvas a ponteado). 

Nota 1: Os sistemas de equações transcenden- 
tes (8), (12') e (16”) resolvem-se com uma certa 
facilidade partindo de pares de valores arbitra- 
dos para x, e 4,. Das primeiras equações extraem- 
-se os valores dos parâmetros (», £ ou 1), nas 
segundas equações determinam-se os respectivos 
i.e.j. completos (K) por meio de tabelas de cor- 
respondências (3, pág. 106/109], e finalmente 
pelas terceiras equações obtêm-se os valores de 
c que correspondem aos pares (xy, 49) inicial- 
mente arbitrados. Como as correspondências 
K (2), K (:) e K (n) são biunívocas, e dado que c 


só pode ser positivo, facilmente se conclui das 
primeiras equações que podem existir três valo- 
res para xy, O que acontece em relação às solu- 
ções em cn para c>1,eemsn para c< 1 (ver 
o gráfico 1). 


Nota 2: Quando ==> 0, x)» 0, o que significa 
que, para sistemas lineares, todo o gráfico se con- 
centra sobre o eixo das abscissas, não perdendo 
em teoria a sua validade quando : =0. É possível 
intensificar as famílias de curvas na região dos 
xo pequenos, ampliando sucessivamente a escala 
das ordenadas até atingir (em teoria com ampli- 
ficação infinita) praticamente a forma das curvas 
de amplificação no caso linear. 


Nota 3: Para cada valor de K (ou do pará- 
metro que lhe corresponde) há uma relação linear 
entre c e x, [ver as terceiras equações de (8'), 
(12º) e (16')]. As famílias de rectas centradas na 
origem que se apresentam no gráfico I, sendo 
linhas de parâmetro constante, são indicativo da 
forma da perturbação (e, por conseguinte, do 
movimento), e para que haja uma noção imediata 
da sua configuração, traçaram-se: três sinusoides 
de parâmetros $=-0, ;- 0,5 e :=-0,9, além da 
curva limite (rectangular) para £==1; e três cosi- 
nusoides de parâmetros «= 0," =0,55e1=-0,9, 
além da curva limite (de impulsos) para u=1. 
Verifica-se que as formas das curvas pouco dife- 
rem das funções circulares até ao parâmetro de 
valor 0,5, esperando-se pequena alteração no 
diagrama na hipótese de perturbação sinusoidal 
pura, na zona compreendida entre as duas rectas 
de «==0,5 e 4==0,5, passando pelo eixo das 
abscissas. 

Também se representou a função dn para 
n==0, n=0,5 e n==0,9 (vindo para n=1 0 
mesmo que cn para « ==1), embora este tipo de 
oscilação não apresente interesse, visto que se 
trata de oscilação provocada por uma força 
de valor médio diferente de zero (a traço inter- 
rompido no diagrama). No caso especial de ser 
n==0 a força será estática, provocando uma de- 
formação também estática do sistema (pontos de 
encontro das curvas a ponteado com as de traço 
interrompido). 


Nota 4: Não abordaremos neste trabalho o 
aspecto da estabilidade das soluções. Este pro- 
blema tem sido satisfatóriamente tratado, (ver 


por exemplo [4, pág. 213 e seg.)). No gráfico 1 
marcaram-se as linhas (a traço interrompido) que 
limitam as zonas de instabilidade (na fig. 5, áreas 
tracejadas). Depende das condições transitórias 
da oscilação, o sistema seguir os ramos estáveis 
superior ou inferior, uma vez que o terceiro 
ramo é instável nas zonas de correspondência 
indeterminada entre c e xo. 

Nas condições reais, uma vez que existe amor- 
tecimento, as curvas fecham-se tal como acon- 
tece nos sistemas lineares, e a amplitude da 
oscilação cai bruscamente do ramo superior para 
o inferior quando se provoca um aumento de 
frequência quando : > 0, ou uma diminuição de 
frequência quando :<0, ou salta bruscamente 
do ramo inferior para o ramo superior se a fre- 
quência desce quando :>>0, ou se sobe quando 
E< 0, assim que se atingem pontos das curvas 
com tangente vertical (fig. 5). 


1 4 


45 
R E >0 
y 
1 c 
RS 
NE Eco 
N 
Fig. 5 


Há que tomar em consideração que, para a 
mesma frequência (tal como estamos a abordar 
este problema), às três soluções xo, correspon- 
dem diferentes valores de K (ou dos parâme- 
tros), e portanto, diferentes formas da função 
perturbadora, porque diferentes serão também 
os coeficientes de elasticidade linear na equa- 
ção (4). 


2.2. Função perturbadora sinusoidal 


A equação diferencial (2) toma a forma 


X 
— + sr x = à cos vt 
Da 


(17) 


q=m=1 


/ Ds : RA 


15 / / 
TT ccencacfecsobrcs lose doa 

CY SR 

A A A 


ARO 


1 E vd 
o A E AO O 


E<O 


deja % | GO | RE 
| ç o 
og E o pd ga ão à % 
da ssiavm trigo: 
| LERU4eELTI 
15 » Ps 
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supondo que a perturbação é puramente sinusoi- 
dal'* de pulsação v==cmw,. Os sinais + ou — 
referem-se respectivamente aos sistemas duros 
(: >> 0) ou macios (:< 0). 

Não se conhece uma solução exacta para esta 
equação. Vários métodos têm sido propostos para 
a sua resolução aproximada. Vejamos resumida- 
mente alguns. 


2.2.1. Método de Duffing **) 


Supõe-se como primeira aproximação que a 
solução é a do sistema linear correspondente 
X == Lo COS CWot (18) 


e substituindo em (17) e tendo em atenção que 


3 1 
cos* cm, t==— cos co, t + nica 3co,t 
4 


obtém-se uma nova equação 


5 sa 
=— (nto 4 ) cos cmo t -- 
4 / 


E 
+ % cos 3 cm;t 


(19) 
que, integrada, permite obter 


1 3 
Ka == =] (o ir q—8,) cos cmo t + 


x, cos 3cm,t (20) 


“36€ 
Nesta primeira aproximação, a família de cur- 
vas xo == f (c,4) será aproximadamente obtida 
comparando os coeficientes dos cos cmot em 
(18) e (20). 


Xy = : (394) É ro 
c* 4cº 


o 3 
4 


ou 


c”) q — do (21) 


ficando (20) com a forma 


= 1 o XÃ 
x2=*,| COS Cont T— pal cos 3 cm, | (22) 
a 36 c 


(*) É indiferente supor a perturbação com uma das 


formas Aosinvt ou A cosvt, não surgindo dai ne- 
nhuma particularização. 


(**) Ver, por exemplo |5, pág. 4-12). 
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As curvas de amplificação serão obtidas para 
cada par de valores (4,c) determinando x, em 
(21) e o valor máximo de xs em (22). 

O método permite maior aproximação. Basta 
reiniciar a substituição em (17), agora com a 
equação (22). 

A função x, = f (c, 4) será agora 


cio Er ) + 
4º 48 


tias (7). 
2x36º À c/. 
e 
1 1 : 
n=n) cos cont + Er rag ( e) edi 
-- 36 49€ c 
1 Xa ' 1 Xn o 
“+ ( el ) eos ses t+ 
924) c 36º À , 
+ x 6 
+| A ( o ) rege qua |cossemt-+ 
25x48 | c 100x<36* c 


es 3 xo Nº 
E cos 7 cont + 
4<49x36º 1 c 


1 za Nº 
Efo e cos 9 copt 
4x 81x 36º" c 


Xo (23) 


(24) 


A iteração poderia continuar indefinidamente. 


Nota 1: Para o mesmo valor de x, a expressão 
(24) difere tanto mais da expressão (22) quanto 
menor for c, quer dizer, quanto menor for a fre- 
quência perturbadora. 


Nota 2: Para a mesma relação de frequências, 
quanto maior for o factor de não linearidade 
(notar que w==V:xo ou V—:x)), maior será 
x, e, portanto, a diferença entre x e Xs. 


Nota 3: Este método iterativo de Duffing só é 


Xy 


aceitável para valores de a suficientemente 


baixos, portanto para sistemas quase lineares 
ou em que a frequência perturbadora é muito 
elevada relativamente à frequência própria do 
sistema linearizado. Acerca da sua validade, 
ver [4]. 


2.2.2. Método de Rauscher (*) 


Este método supõe como primeira aproxima- 
ção que, para uma determinada amplitude de 


(*) Ver, por exemplo, [5, pág. 4-13). 
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oscilação, a frequência é a natural do sistema 
livre correspondente. 
No caso presente 


x 


=+rtxté=o0 
Do” 


admite soluções com a forma [1]: 
— para :>0 (sinal +) 


x=xmcn(u|p) 


com 
p? 
pesa, — a 
2 + 2x5 (26) 
e 
aca / 2 
u=9)V1+4 % t 
— para *< O (sinal —) 
x=xsn(ulé) 
com 
2 — x) 
e 


Os períodos de oscilação correspondentes serão 


— para E > 0 
E: 4K (u) 
Am Va (27") 
AM > 
— para E< 0 
OE. K (2) 
po 
n Va (27º) 


1 
= 98 
2 


Para que o sistema oscile com a frequência 
imposta »j = ci w, terá que ser 

2% 

VI cm Ci g = ——— 


T (28) 


Através das expressões (26) e 28) obtemos em 
seguida a função t (x) 


— para * >0 


t == 


x 
; — aTg Cn —— 
o) | j «+ ão Xy 


(29º) 


— para E (O 


vindo 
A, cos co t = À, cos e +) | 
T 
e a equação (17) fica 


É +(2— a, cos[ 25102 |) + 4 =0 
OT) Ti 
(30) 


que é uma equação autónoma, como (25) tam- 
bém da forma «x + f(x) =0, a qual admite uma 
solução t (x) de período Ts: e frequência angular 


2” 
va==C20, = —— 


E (31) 


Substituindo novamente em (17), obtinha-se 
assim um método iterativo. Este método é mais 
rapidamente convergente que o de Duffing, mas 
tem o inconveniente de originar equações dificil- 
mente integráveis. No caso em estudo, a equa- 
ção (30) só poderá ser integrada por métodos 
gráficos, ou das diferenças finitas. 

Notar que no método anterior nos fixámos 
num par de valores (4, e c) e fomos sucessiva- 
mente corrigindo o valor de xo. Neste fixou-se 
A,e x,, e vai-se corrigindo c por iteração. 


+ + 
Baseado neste método proporemos uma va- 
riante '*) que nos parece ser mais aproximada. 
Em vez de considerar a equação (25) começa-se 


por usar a que se apresentou em (4) vindo para 
soluções 


— para E >0 


Z = go cn (u/p) 


(*) Esta variante tanto se aplica a sistemas deste tipo 
como a outros quaisquer não dissipativos pertubados por 
quaisquer funções p (t), vindo (32) na generalidade 


x+[02x—p(t6)]+Ht(G)=0 


em que t (x) se obtém integrando a equação especial do 
movimento com perturbação de forma favorável, 


n , o x (f(x) não contém 
x + ogx—+ E (x) =p (t) = po — termo linear) 
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com 
x 
p=—— (32) 
24,—-28,+4+ 2% 
AR 
u = v04/ 1—— Va + mx t 
o 
e 
T— AKI) 
/. A. Í 9 
Do V o MV 1 + No 
— para E< 0 
q = a, sn (u/t) 
com 
sp 
, = | (327) 
2%—28,—4, 
a / 1 «2 
U= 6% V a SR V O Eta: t 
e 
4 K (i 
= —— E fo 


Tendo em conta (28), a função t(x) virá 


— para : >0 
es a x 
= 7 arg cn x (83) 


Do V 1 — j V à 1 + Xo 


+ T/4 “0 
T=4 | dt=4 


Ti x 
a Ses O 
4K (g) No 
— para E< 0 
1 x 
= ———————— o e arg sn ——- = 
O Va — da 1 aca de is 
ha Xo 24 
(337) 


a 
Ko SO a 


Segue-se tal como anteriormente a expres- 
são (30), mas agora com t (x) obtido já bas- 
tante próximo do resultado exacto (fig. 6). 

A iteração continua, encontrando-se as mes- 
mas dificuldades na integração que no método 
de Rauscher. 


Estas operações são bastante trabalhosas e a 
sua resolução num calculador automático obriga 


a uma programação bastante complexa devido à 
necessidade de memorizar as tabelas de funções 


elípticas ou as séries correspondentes e = 
Xa 
x x 

=snu oucnu,eu=argsn-—— ou arg en — ). 
x Y 
nad xy 


Manualmente, resolveram-se 4 pontos com 
uma iteração, através da integração numérica da 
expressão 


0 “ Xo 


Yo 


ao 


A 


ME TODO DE RAUSCHER 


[, q 1 . Pa T E 
nós) (mw -a')p 24, [) Cos] Trarg en  | da 


procurada 


x, Sin | 2K sn X 


MÉ TODO PROPOSTO 


Fig. 6 
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obtida pela substituição de (33) em (30), e 
integrando estas duas vezes. Os sinais + ou — 
e as funções cos ou sin, e cn ou sn, usam-se 
quando é :>>0 ou :<[0, respectivamente. O novo 
valor de c virá dado por 


sultados obtidos são pouco aproximados (ver 
quadro 1). Os valores corrigidos estão assinala- 
dos no gráfico 1. 


QUADRO I 


27% - 
Cc = To pes; aja spa Solução corrigida 
"0 e | ce 
do X c Ar | Na C 
Consideram-se intervalos 4x == 0,1, tomando- | — 
-se valores médios na resolução numérica do :>0 1,5 |1,195| 0,88| 1,5 asá 0,84 
primeiro integral. Na integração final efectuou-se so) —1 li | 1,67] —-1 l111 | 1,69 
a resolução gráfica, uma vez que tende e<O0| 0,2 |0,78 | 0,50: 0,2 |0,78 | 0,54 
|<o —7 11,055 | 1,04 — 71 |1,055 1,06] 
assintoticamente para co quando x tende para | 
zero. Além disso, as operações foram efectuadas 
com régua de cálculo, e por consequência os re- (Continua) 
SUMÁRIO 
Apresentam-se alguns métodos para o estudo das amplitudes de oscilação de sis- 
temas com um grau de liberdade, não amortecidos, do tipo Duffing, perturbados por 
forças sinusoidais circulares e elípticas. 
SYNOPSIS 
The author presents some methods for studying vibration amplitudes of undam- 
ped one-degree-of-freedom systems, of the Duffing type, acted by circular and ellipti- 
cal sinusoidal perturbatious. 
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NOTAS INFORMATIVAS C. D. U. 621.341.5/94 
Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 
— Elementos extraídos das estatísticas mensais do Repartidor Nacional de Cargas (R. N. C.) 


Nora: As produções e os consumos das empresas do R, N. C. representam 
I — Breve nota mensal cerca de 94 º/, dos totais do Pais, 


No aspecto hidrológico, o mês de Junho apresen- 
tou-se, no conjunto, muito acima da média. 


Consumo em bombagem (Chb)(!) 
Produção para con- 
sumo não perman, (Penp) (2) 


II — Elementos gerais (GWh) JUNHO 
a) Mensais 
riação 
1963 [ooo | ass | Mo 1964 | Jatine | 
lo 
Produção hidráulica Vi 313,2 | 352,8 | +- 13 AR 
Produção térmica (P+)..... 0,0 0,0 O ias 
Produção total (PT). ..... 318,2 | 852,8 | 18 A 
Exportações (Ex)... ...... 0,0| 0,0 Amis nua e 
Importações (1) . ........ 0,0 0,0 amis À aco HH 
Consumo em bombagem (Cb)(!)) — — — na ARNnam 
Produção para con- | REDE RAGEE 
sumos não perman. (Penp)(2) | 16,0) 81L0/+ 6 nunadránam 
Produção para con- DERA Todd 
sumos permanentes (Pep) (3) 236,1 211, 8 | +14, 8 (1 SESECÓREA 
Coeficiente de hidraulicidade| 1.20) 1,25 — tt 
fa Ei GG E fã 
b) Acumulados desde 1 de Janeiro EESEGRSEE 
Variação HEEHEHHHA 
1963 | 1964 0/ " ú tá Ss - 
U 2a 
EE mes ETR —) [ur Es] o 
Produção hidráulica (Ph)... 2049,9 /2236,5| + 9 
Produção térmica (P+)...... 18,9) 594/4214 
Produção total (PT)..... ek al qa + e 
) DR st ace emo 32 36,4 i 
Rb th ento ' 0.1 0,0 a 100 IV — Energia armazenada nas principais albufeiras 


438,4 | 478,9 |. 9 No fim do mês 


Albufeiras : ES 


Produção para con- | 
sumos permanentes (Pcp) (?) 1597,7 |1781,8 | 4-11,5(1) GWh | 0, (1) 
Coeficienta de hidraulicidade | ] Yi 1,11 | cu 
| 
NOTAS 02 ( é 
(1) Energia consumida em bombagem para armazenamento na Al- Alto Rabagão o e MR DO NRO» i 194, 20,5 
bufeira do Alto Rabagão. aaa FAS das Es Es 201,4 40,9 
2) Conside - u ã tes os indic 2. RES + o 
(3) A rodGas Eerk cónenaD pemianenias (bc) É otra Venda Nova « «cerva. 94,5 141 
pela seguinte expressão: Pcrp —- P7—-Poop/ll—Ex Salâmondo . .« cc cv vc vo 21,0 76,1 
(4) O aumento percentual de produção para consumos permanen- Caniçada o: dt & e mude: E 31,4 94,5 
tes, tendo em conta a incidência q dos domingos e dias especiais, é Vilar 
respectivamente de 10,8º/, e 10,9º : ai RO De A O dA “om ] peca — 
IL — Diagramas de carga tá dias característicos Cabrili . scr. 668,8 19,2 
- jon Castelo do Bode. . .,..| 1441 88,4 
4.º feira; Gullhonel «css cc Ra 7,9 95,2 
1963 19614 Lagoa Comprida +... ..| 24,9(?) 73,5 
per" ti AO o ado sou só 09,5 
Produção hidráulica (Ph) MWh/ 11390 | 12380 ii ia Po 67,8 
Produção térmica (Pr) MWh O 0 FERA FS Ea da: NU , 69,0 
Produção total (Pr) MWh| 11390 | 12380 Póvoa. cc wcca ros 9,4 (*) 69,2 
Trocas com [ Export. (1x) MWh E j Total | COM A. Rubagão . « «| 1065,7 — 
Espanha Import. (1) MWh abas O ' 
TOTAL PT + (L-Ex) MWh | 11390 19380 sem A, Rabagão . . .| 866,2 83,0 
Prod. para cons. perm. — (Pcp) MWh 9043 J146 
Prod. para cons, não perm, (P ie) Mwh| 2947 | 2634 + 
Potência máx. MW 626 bTo É 
2e|p. LEx) Potência mín. MW 310 3359 Nodas: 
SE THEO  Utilizda ponta horas | 18,2 183 1 4 
% 4 Factor de carga 0,76 0,76] (!) Coeficiente de enchimento em energia das albufeiras, 
sz ia o == = —uw| 520 567 | (2) Inclui 3,3 GWh armazenados em Vale do Rossim no 
É - A Potência min. MW 208 224 início do mês e 3,3 GWh no fim do mês, 
Us ep Utiliz. da ponta horas AR ca (3) Inclui 2,5 (GWh armazenados no açude do Poio no 
| Factor de carga 0,12 0,f início do mês e 2,2 GWh no fim do mês, 
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RESERVATÓRIOS CILÍNDRICOS DE BETÃO 
ARMADO PARA ÁGUA 


CÁLCULO GRÁFICO E PROPORÇÕES ECONÓMICAS 


1 — INTRODUÇÃO 


O presente trabalho, embora não devesse ser, 
foi elaborado por um simples particular sem re- 
muneração e sem possibilidades de recurso aos 
computadores electrónicos e aos aparelhos de 
preparação de gráficos. 

Espero, com a sua apresentação, estimular o 
aparecimento de outro melhor que o substitua e 
complete, elaborado pelos competentes Serviços 
Oficiais. 

O trabalho é dividido em duas partes a saber: 

Cálculo gráfico 

Proporções económicas . 


2 — CÁLCULO GRÁFICO 
2.1 — Generalidades 


Revistos os métodos de cálculo de diversos 
autores, a atenção concentrou-sesobre umtrabalho 
de M. Hangan e M. Soare intitulado «Calcul Ra- 
pide des Reservoirs Cylindriques» (1) embora 
nele se trate apenas o caso de reservatórios com 
parede livre em cima e soleira inextensível e 
assente em fundação rígida, sendo a espessura da 
parede igual em toda a altura e a da soleira também 
constante para qualquer ponto da superficie. 
Os autores apresentam as fórmulas que conduzem 
à solução do problema e uma série de ábacos nos 
quais se entra mediante o cálculo dos seguintes 
parâmetros: 


A BH='*y3(1—4) 15 
VRe 

2 ee 

em que: 


R ==raio da parede cilíndrica 

H== altura da água e da parede cilíndrica 
e =espessura da parede 

» soleira 


Í 
É = » 
v == coeficiente de Poisson 


ANTÔNIO MANUEL DE VASCONCELOS 


Eng.º Civil (L.S.T.) 


Com estes dois parâmetros os ábacos dão uns 
coeficientes que multiplicados por wH', vRH e H 
(sendo » o peso específico do líquido) dão os 
momentos flectores, as tracções e as cotas dos 
pontos especiais. 

Por considerar este caminho demasiado longo 
construi novos ábacos nos quais se entra directa- 
mente com os dados do problema e um produto 
fácil de calcular — Re — e onde se lêem imediata- 
mente os resultados sem necessidade de qualquer 
cálculo além do já citado produto Re. O agrupa- 
mento destes ábacos, como se apresenta, permite 
mesmo eliminar a necessidade de se fazer refe- 
rência às cotas de alguns dos pontos especiais 
e permite ainda construir sobre os próprios 
ábacos os diagramas das tracções e dos momentos. 

Além disto procurou-se resolver o problema de 
forma a considerar o caso mais frequente da 
espessura da soleira na proximidade da parede 
ser maior do que na zona central. 

Como este caso de soleira de espessura variável 
não se encontra tratado teoricamente no estudo 
citado, procedeu-se da seguinte maneira: 


— No cálculo das tracções e momentos na 
parede (excepto na ligação à soleira) admite-se 
que a espessura da soleira é constante e igual à 
da parede porque é o caso mais desfavorável e 
que, portanto, dá maior segurança. 

— No cálculo do momento de encastramento e 
da tracção na soleira, admite-se que a espessura 
da soleira é também constante mas dupla da es- 
pessura da parede, caso mais desfavorável do 
que se as espessuras fossem iguais conduzindo, 
portanto, a maior segurança. 


Na construção dos novos ábacos tentou-se 
considerar o alongamento radial da soleira, pois é 
o caso real e o que conduz às maiores tracções e 
aos maiores momentos flectores na parte central 
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da parede. Porém, as fórmulas apresentadas no es- 
tudo citado não são suficientes para este fim, 
pois, apenas é dada a fórmula para o cálculo do 
momento na base da parede, que é exactamente 
onde menos interessa, porque aí, o momento 
considerando a soleira extensível é menor do que 
se esta for considerada inextensível. É certo que 
a diferença não é grande (7!/y no exemplo apre- 
sentado naquele estudo) mas ignoro se não será 
importante nas tracções e nos momentos da zona 
central da parede. 

Assim, os novos ábacos, como os do estudo 
citado, são estabelecidos na suposição irreal de 
que a soleira não se alonga radialmente. 

Dada a enorme soma de trabalho que implica- 
ria o estabelecimento dos novos ábacos a partir 


a)-TRACÇÕES 


MN No No NM 


Diagrama exacto 


-—— — — — Tinha auxiliar 


aÃ Diagrama adoptado 


das fórmulas, algumas bastante complexas, uti- 
lizaram-se valores lidos cuidadosamente nos ába- 
cos do estudo citado. 

Houve, no entanto, que fazer os cálculos para 
os valores correspondentes a LH=-4 que está 
fora do âmbito destes ábacos mas que me pare- 
ceu interessante considerar. 


2.2 — Tracção nos anéis 


Como já foi dito a tracção nos anéis é consi- 
derada supondo a espessura da soleira constante 
em toda a superfície e igual à da parede. 

O ábaco preparado (n.º 1) é constituído por 
três ábacos acoplados e tem adjacente um es- 
paço para o diagrama dos esforços. 
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Fig. 1 — Diagrama das tracções e dos momentos nas paredes livres em cima e ligadas às soleiras 
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ÁBACO N.º 1 


Fig. 2 — Tracção nos anéis duma parede livre em cima e ligada a uma soleira de fundação rígida 


NOTA - Os valores de N referem-se a Iml. de altura de parede 


Mediante o conhecimento de R, de He dee 
bem como do produto Re, obtém-se os valores 
da tracção unitária nos diferente anéis em que se 
considere dividida a parede, sem necessidade de 
qualquer cálculo intermédio. 

Por não ser prático traçar o diagrama exacto 
dos esforços, propõe-se a sua substituição por 
duas linhas auxiliares conforme se indica na 
fig. Ia passando cada uma pela sua extremidade 
da parede e ambas pelo ponto de tracção má- 
xima dado pelo valor de c (do ábaco superior 
direito) e de N (do ábaco inferior esquerdo). 

Considerando-se suficiente uma análise por 
faixas de 0,5 m de altura, o campo do diagrama 
foi marcado com os respectivos traços pelo que 
não será necessário desenhar as linhas auxiliares 
acima referidas, pois, uma régua ou esquadro 
colocado nas suas posições, permitirá ler os va- 
lores de N; correspondentes a cada faixa fi; e, 
assim, adoptar um diagrama em degraus con- 
forme se indica na fig. la. A faixa superior fi- 
cará com uma altura diferente de 0,5 mse H 
não for múltiplo de 0,5 m. 

Para não se esquecer que os valores de N se 
referem a 1 m.1. de altura da parede, apresenta-se 
uma sugestão para o quadro de valores a ela- 
borar : 


N.º de Atica Tracção 
ordem unit. 
sã da É 
faixa max. na 
faixa fi faixa 
i Ni 
tm) | (kg) 


Na escolha das armaduras deve ser-se gene- 
roso, pois, é necessário não esquecer que o ábaco 
foi organizado na suposição de que a soleira era 
inextensível, quando, na realidade o não é e, por- 
tanto, as tracções são superiores às indicadas. Em 
contrapartida, porém, a maior espessura da so- 
leira junto à parede conduziria a menores trac- 
ções. 

Embora tenha sido preparado para reservató- 
rios de água, de peso específico » = 1000 kg.m?, 
ele pode ser usado para qualquer outro líquido 
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(Ra=: 


desde que se multipliquem os valores de N pela 
densidade deste líquido. 


2.3 — Momentos flectores na parede 


Como já foi dito apenas se considera o caso da 
soleira ter espessura constante e igual à da parede. 

O ábaco preparado (n.º 2) é constituído por 
dois ábacos acoplados e tem adjacente um espaço 
para o diagrama dos momentos. 

A utilização é analoga à do ábaco n.º 1, apenas 
havendo que fazer uma referência à forma do 
diagrama dos momentos. O traçado da linha 
auxiliar da zona inferior da parede não convirá 
fazer-se por AB mas sim por AC (fig. 1b). Com 
efeito, quando uv reservatório estiver vazio e o 
terreno não for perfeitamente rígido, haverá trac- 
ções na face exterior da base da parede e na face 
inferior da soleira cujos valores, variarão conforme 
as características do terreno. Além disto as arma- 
duras verticais da zona central da parede convém 
serem prolongadas até à base para facilitar man- 
ter em posição as armaduras anelares durante a 
construção. 

No que respeita à divisão da parede em faixas 
de 0,5m de altura, à generosidade na escolha 
das armaduras e à possibilidade de utilização 


Armadura adoptada 


Secção pi 
da Nºeg Secção 
armadura dos efectiva 
necessária varões 
Nixf 
Aai= rs 
Ra 
« « kg/em?) | 
(cm?) (mm) | (cm?) 


para outros líquidos que não a àgua, tem aqui 
lugar as mesmas observações que se fizeram 
para as tracções nos anéis. 

Embora a compressão proveniente do peso da 
parede e da cobertura e suas cargas seja ge- 
ralmente pequena, não se deve esquecer o seu 
efeito nas tensões da flexão. 

Na maioria dos casos adopta-se em toda a 
altura da parede a mesma armadura, mas, no 
caso de grandes reservatórios, pode convir dife- 
renciar a armadura por faixas de altura e, então 
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ÁBACO N.º 2 


Fig. 3 — Momentos flectores numa parede livre em cima e ligada a uma soleira de fundação rígida 


poderá preparar-se um quadro com a seguinte 
disposição: 


varia com H, conforme se pode observar pelo 
seguinte quadro, com valores relativos ao caso 


SS 


Armadura 
N.º de Mom. M; Compressão no betão 100 Aaji |——— À a 
ordem flect. qe R'bi b'h Secção | Varões 
da bh? Ra = ke./ > necess. - 
Mi lh =...) (Ra =.. .. kg. /em*) Aai N.' Secção 
faixa (hê = ) i K j efectiva 
(kg:m) E (kgs /em?) (em?) | (mm) | (em?) 


O ábaco dos valores de b (da direita) resultou 
pouco expressivo porque a variação de b é muito 
pequena em relação a H. Poderia ter-se, até, 
suposto b constante para qualquer valor de H, 
análogamente ao que se fez nos ábacos n.º* 3, 
4e5. 


2.4 — Momento flector na ligação parede-soleira 


Como se disse o caso mais vulgar será o da 
soleira ter espessura dupla da parede pelo me- 
nos na base desta. Apresentam-se, no entanto, 
os ábacos referentes aos três casos seguintes: 


ee' =1 ÁAbaco n.º 3 
ele == 0,5 » » 4 
ee'=0 » » O 


O primeiro caso pode interessar para peque- 
nos reservatórios em que não se preveja o re- 
forço do ângulo parede-soleira. 

O segundo caso corresponde a uma soleira 
de espessura uniforme dupla da espessura da 
parede mas deve conduzir a valores pouco supe- 
riores aos do caso desta sobrespessura da soleira 
se verificar apenas numa faixa junto à parede, 
que corresponderá à forma mais vulgar. 

O terceiro caso corresponde ao encastramento 
perfeito e é apresentado apenas a título de 
curiosidade para mostrar os valores limites, pois, 
na prática, a parede nunca será perfeitamente 
encastrada na soleira. 

Qualquer destes três ábacos é constituído por 
três ábacos acoplados. Os da esquerda dão os 
momentos em função de H e de Ke. Os da direita 
dão a altura da parede designada por a, em que 
o momento se anula supondo um valor de 
H= 55 m. Procedeu-se desta forma porque 
trouxe simplificação dos desenhos e porque não 
tem inconveniente. Com efeito o valor de a pouco 
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de encastramento perfeito, que é o mais desfavo- 
rável: 


Re H a 
(m?) (m) (m) 
0,1 3,6 0,18 

» 5,5 0,21 

1,4 3,6 0,58 

' 5,5 0,64 


Os ábacos inferiores dão a distância em que 
a soleira é afectada pelo momento. Usou-se a 


fórmula d=21/ E devida a K. Girkmann (2) 


em que P é a carga da soleira. Os valores de P 
foram calculados para água e para um peso da 
soleira variável conforme a sua espessura prová- 
vel. Análogamente ao referido para o ábaco an- 
terior, supõe-se d constante para qualquer valor 
de H, donde não vem inconveniente de maior, 
conforme se pode observar do seguinte quadro, 
com valores relativos ao caso de encastramento 
perfeito, que é o mais desfavorável. 


Re H d 
(m”) (m) (m) 
0,1 3,0 0,31 

» 5,5 0,32 

Lo 3,0 1,05 

» 9,9 1,14 


Ao contrário do que acontece com as tracções 
nos anéis e os momentos flectores na zona cen- 
tral da parede, não se julga necessário traçar os 
diagramas dos momentos na zona da ligação 
soleira-parede porque a extensão das armaduras 
é geralmente pequena. Julga-se, no entanto, de 
sugerir que as armaduras sejam reduzidas a me- 
tade na metade superior da faixa de parede de 
altura a e na metade interior da coroa circular 
de largura d. Nesta última zona haverá mesmo 
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